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RESPONSE  OF  AN  ACOUSTIC  COATING  ON  A  RIB-STIFFENED  PLATE 


1.  INTRODUCTION 


Underwater  vehicle  hulls  are  typically  metallic  with  acoustic  coatings  that  are  made  of  a 
polymer  material.  These  hulls  are  usually  internally  reinforced  in  one  direction  to  provide 
increased  stiffness  against  the  hydrostatic  forces  that  act  on  the  hull  when  the  vessel  is 
submerged.  Acoustic  coatings  are  used  to  produce  quieter  vehicles,  contain  sonar  sensors, 
inhibit  drag,  and  prevent  biological  fouling  on  the  exterior  of  the  vehicle.  Understanding  the 
structural  response  of  such  a  system  is  important  for  the  design  of  new  underwater  vehicles  and 
analysis  of  existing  underwater  vehicles.  An  analytical  model  of  this  type  of  configuration  will 
allow  predictions  of  the  dynamic  response  of  an  acoustic  coating  on  a  marine  structure  and  any 
sonar  system  that  is  embedded  in  the  coating. 

The  research  of  plate  theory  has  been  an  ongoing  field  of  study  for  many  years.  Thin  beam 
(Euler-Bemoulli)  theory  originates  in  the  eighteenth  century.1  This  model  is  sometimes  called  a 
flexural  wave  model,  and  an  extremely  similar  model  exists  in  plate  theory.  Mindlin2  modified 
thin  plate  theory  to  include  the  dynamics  of  rotary  inertia  and  shear  effects,  and  this  extended  the 
low  wavenumber  accuracy  of  the  model.  Analysis  of  wave  propagation  in  fully-elastic  plates 
consisting  of  one  or  multiple  layers  using  Naviers’  equations  of  motion  has  been  studied 
extensively  and  is  well  documented.2  Fluid  loading  has  been  added  to  thin  plate  analysis4  and 
thick  plate  single  and  multilayer  analysis.5  Using  these  techniques,  the  dynamic  response  of 
fluid-loaded  acoustic  coating  on  an  unreinforced  plate  can  be  calculated  and  analyzed. 

Over  the  years,  plate  theory  has  been  expanded  to  include  stiffening  effects  of  ribs  that  are 
typical  of  many  marine  and  aviation  structures.  They  are  almost  inclusively  placed  on  one  side 
of  the  plate  to  reinforce  the  structure.  The  response  of  a  periodically  supported  beam  to  a  load 
with  fixed  wavenumber  and  frequency  has  been  investigated.6  The  response  of  periodically 
stiffened  fluid-loaded  plates  to  harmonic  loading7  and  line  and  point  forces*  has  been  established. 
The  problem  of  aperiodicity  in  the  stiffeners  has  been  solved.4  This  problem  was  also 
investigated  for  a  finite  number  of  equally-spaced  stiffeners10  and  randomly  spaced  stiffeners.11 
Asymptotic  models  of  plate  radiation  into  fluid  fields  have  also  been  developed.12  It  is  noted  that 


these  papers*'  12  use  some  form  of  thin  plate  theory,  and  the  resulting  frequency  limit  where  the 
model  assumptions  are  valid  depends  on  the  thickness  of  the  plate,  but  for  most  plates  it  is 
typically  in  the  hundreds  of  Hertz.  The  acoustic  response  of  a  fully-elastic  cylindrical  shell  with 
a  complete  acoustic  coating  has  been  researched12  as  well  as  that  of  a  cylindrical  shell  with  a 
partial  acoustic  coating.14  There  are  no  rib  effects  present  in  either  of  these  papers.  Higher 
frequency  analysis  is  possible  using  numerical  methods  such  as  finite  element  analysis;15 
however,  these  computations  can  be  time  consuming  and  frequently  have  stability  problems 
when  Poisson’s  ratio  of  the  coating  material  approaches  0.5.  Recently,  a  fully-elastic  solution  to 
a  plate  containing  discrete  masses  was  developed.1*’  This  method  can  be  extended  to  model  the 
behavior  of  an  acoustic  coating  on  a  ribbed  backing  plate. 

This  report  presents  an  analytical  model  of  a  plate  in  contact  with  a  fluid-loaded  acoustic 
coating  on  one  side  and  a  series  of  equally  spaced  ribs  on  the  other  side.  The  structure  is  loaded 
through  the  fluid  with  an  incoming  acoustic  wave.  The  plate  and  the  acoustic  coating  are 
modeled  as  three-dimensional  fully-elastic  solid  bodies,  the  fluid  is  modeled  as  a  three- 
dimensional  acoustic  field,  and  the  ribs  are  modeled  as  discrete  springs.  The  formulation  begins 
with  elasticity  theory  where  the  motion  in  the  plate  and  the  acoustic  coating  are  modeled  as  a 
combination  of  dilatational  and  shear  waves.  These  waveforms  can  be  used  to  determine  the 
three-dimensional  displacement  fields  with  unknown  coefficients.  These  displacements  are 
inserted  into  stress  and  continuity  equations  at  the  system  boundaries  and  interfaces  that  contain 
the  system  excitation,  the  fluid  loading,  and  the  force  of  the  ribs  on  the  structure.  Using  an 
orthogonalization  procedure  produces  an  m-indexed  mathematical  model  of  the  system  where 
each  m  index  is  a  set  of  equations  twelve  rows  by  an  infinite  number  of  columns.  All  of  the  m- 
indexed  equations  can  be  combined  and  this  yields  a  matrix  system  of  infinite  extent,  which  is 
truncated  to  a  finite  number  of  lower  wavenumber  terms.  Inverting  this  matrix  solves  for  the 
unknown  coefficients  and  produces  system  displacements,  stresses,  and  pressure  field  solutions. 
The  fully-elastic  model  is  compared  to  a  previously  developed  thin  plate  model  at  low  frequency 
and  wavenumbers  to  insure  accuracy  and  consistency  to  prior  work.  A  numerical  example 
problem  is  then  studied  with  specific  interest  in  reception  of  acoustic  signals  and  strength  of  the 
scattered  field.  It  is  shown  that  the  ribs  can  moderately  affect  the  system  response. 
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2.  SYSTEM  MODEL 


The  system  model  is  that  of  a  fluid-loaded  solid  layer,  called  the  acoustic  coating,  in 
contact  with  a  rib-stiffened  solid  layer,  called  the  backing  plate,  as  shown  in  figure  1 .  This 
problem  is  analytically  modeled  by  assuming  that  the  pressure  in  the  fluid  is  governed  by  the 
acoustic  wave  equation,  both  solid  layers  are  governed  by  fully-elastic  equations  of  motion,  and 
the  rib  stiffeners  are  modeled  as  discrete  linear  springs.  The  acoustic  coating  is  loaded  on  the  top 
surface  with  an  incoming  acoustic  wave.  The  ribs  on  the  bottom  of  the  backing  plate  are  equally 
spaced  at  a  distance  of  L  (m)  in  the  x-direction  and  each  has  stiffness  per  unit  length  of  K  (Nm:). 
The  model  uses  the  following  assumptions:  (1)  the  forcing  function  acting  on  the  plate  is  a  plane 
wave  at  definite  wavenumbers  in  the  x-  and  y-directions  and  frequency  in  time,  (2)  motion  in 
both  plates  is  normal  in  the  z-direction  and  tangential  in  the  x-  and  y-directions  (three- 
dimensional  system),  (3)  both  plates  have  infinite  spatial  extent  in  the  x-  and  y-directions,  (4)  at 
the  interface  of  the  acoustic  coating  and  backing  plate,  the  surfaces  move  in  contact  with  one 
another  and  the  stresses  are  equivalent,  (5)  the  ribs  are  modeled  as  springs  that  have  translational 
degrees  of  freedom  in  the  z-direction,  (6)  the  particle  motion  is  linear,  and  (7)  the  fluid  medium 
is  lossless. 


ZA 


Figure  1.  Schematic  of  Plate  System  Showing  Coordinate  System 
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The  motions  of  both  the  backing  plate  and  the  acoustic  coating  are  governed  the  Naviers’ 
equations  of  motion  written  in  vector  form  as 


//V 2 u (x,y, z,t)+(A  +  //)■ V V  .  u(*, y,  z,t)=pd^UiX'^ZJ) 

dr 


(1) 


where  p  is  the  density  (kg  m  3),  A  and  p  are  the  complex  Lame  constants  (N  m  :),  /  is  time 
(seconds),  •  denotes  a  vector  dot  product,  and  u(.y,v^,/)  is  the  three-dimensional  Cartesian 
coordinate  displacement  vector  and  is  written  as 


u  (x,y,z,t)  =  [u(x,y,z,t)  v(x,y,z,t )  w(x,y,z,t)]' 


=  V</>(x,y,z,t)+S7  xip(x,y,z,t) , 


(2) 


where  ^is  a  dilatational  scalar  potential,  V  is  the  gradient  operator,  x  denotes  a  vector 
crossproduct,  and  ip  is  an  equivoluminal  vector  potential  given  by 


ft(x,y. 


t)=[ipx(x,y,z,t)  ip  v(x,y,z,t) 


ip:(x,y,z,t)]T . 


(3) 


Additionally,  the  constraint  equation  on  the  equivoluminal  vector  potential  is 
V  •  ip(x,y,z,t)  =  0. 

Expanding  equation  (2)  gives  the  individual  displacements  as 

d<f>(x,y,z,t)  dipz(x,y,z,t)  dip  v(x,y,z,l) 
u(x,y,z,t)  = - r - +  — =  ~ - — ; - , 

ox  cy  oz 

..  _  ^_d</>(x,y,z,t)  dipz(x,y,z,t)  dipx(x,y,z,t) 

v\x,y,z,t)  - - - - - - + - - - , 

cy  dx  dz 

and 


(4) 


(5) 

(6) 
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(7) 


,  ^  _d</>(x,y,z,t)  dif/y(x,y,z,t)  di//x(x,y,z,t) 

wyx^y9zj)  - - - - r ■ 


dz 


dx 


dy 


Equations  (5),  (6),  and  (7)  are  next  inserted  into  equation  (1),  which  results  in  four  decoupled 
wave  equations,  given  by 


2^2,,  #v  dz</>{x,y,z,t) 
cdS7  </>(x,y,z,t)  = - - - 

dt 


(8) 


dt2 


(9) 


2  2  ,  d  V y(x,y,z,t ) 

cjV  y/v(x,y,z,t)  = - ^ — - 

dt2 


(10) 


c,2vV;(v.V)-8  r-(Y,,',) 

dt 2 


(ii) 


where  equation  (8)  corresponds  to  the  dilatational  component  and  equations  (9),  (10),  and  (11) 
correspond  to  the  shear  component  of  the  displacement  field.  Correspondingly,  the  constants  cj 
and  cs  are  the  complex  dilatational  and  shear  wave  speeds  (m  s'1 ),  respectively,  and  are 
determined  by 


cd 


X  +  2/j 
P 


and 


(12) 
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(13) 


The  system  equations  of  motion  are  formulated  as  a  boundary  value  problem  using  twelve 
equations  of  stress  and  displacement  continuity  written  in  terms  of  the  plates'  displacements  and 
corresponding  forcing  functions.  The  system  is  loaded  on  the  bottom  by  the  springs  acting  in  the 
z-direction;  thus,  the  normal  stress  at  this  location  (z  =  a)  is  written  using  a  force  balance 
between  the  springs  and  the  lower  plate  as 


rzz(x,y,a,t)  =  A  i 


dux(x,y%a.t)  |  dvx(x,y,a,t)  dw^y^aj) 
dx  8y  dz 


+  2/q 


du’i(.v,y,  aj ) 


n  =  +  oo 

=  ^Kwifay^rfSix-nL)  , 

n  —  —oo 


(14) 


where  S(x-nL)  is  the  spatial  Dirac  delta  function  that  distributes  the  force  of  the  springs 
discretely  and  equally  spaced  in  the  x-direction  and  the  subscript  1  corresponds  to  the  backing 
plate.  The  tangential  stresses  at  the  bottom  of  the  backing  plate  are  zero  and  are  written  as 


Tzv(x,y,a,t )  =  H\ 


dv\  (-TV,  a ,  0  dw,  (x,  y,  a,t) 


dz 


dy 


=  0 


(15) 


and 


rzx{x,y,a,t)  =  v  | 


dw\(x,y,a,t)  |  dux(x,y,a,t) 
dx  dz 


=  0. 


(16) 


The  interface  between  the  backing  plate  and  the  acoustic  coating  is  modeled  using 
continuity  of  displacements  and  stresses.  Continuity  of  the  stress  fields  at  the  interface  (z  =  b) 
yields 
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(17) 


du\(x,y,b,t)  t  dv^yyyAt)  |  y,b,t) 

dx  8y  dz 

dujiyyyJxj)  |  dv2(;r,y,/>,0  |  dwgj^yAt) 
dx  dy  dz 


+  2/u\ 


+  2/j2 


dw^tyyjyt) 

dz 

dw2(x,y,b,t) 


dz 


M\ 


dvx(x,y,bj) 

,  dw}(x,y,b,t) 

dv2(x,y,b,t ) 

|  dw2(x,y,b,t) 

dz 

ay 

=  Ml 

dz 

dy 

and 


M\ 


dwii^yjyt)  |  dux(x,y,b,t) 
dx  dz 


=  Ml 


d^y^yAt)  |  di^yyyAt) 
dx  dz 


(18) 


(19) 


The  subscript  2  corresponds  to  the  acoustic  coating.  Continuity  of  the  displacement  fields  at  the 
interface  yields 


ux(x,y,b,t)  =  u2(x,y,b,t). 

(20) 

v\(x,y,b,t)  =  v2(x,y,b,t) , 

(21) 

w\  (x,  y,b,t)  =  w2(x,  y,  b,  t) . 

(22) 

The  top  of  the  acoustic  coating  is  in  contact  with  a  fluid,  and  the  normal  stress  at  this 
location  (z  =  c)  is  written  using  a  force  balance  between  the  pressure  in  the  fluid  and  the  coating 
as 


rzz(x,y,c,t)  =  A 2 


du2(x,y,c,t )  |  gv2(x,y,c,/)  |  dw2(x,y,c,t ) 


dx 

=  -Pa(x>y,C’t)  , 


dy 


dz 


+  2ju2 


dw2(x,y,c,/) 

dz 


(23) 
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where  pa(x,y,c,t)  is  the  pressure  field  in  contact  with  the  top  of  the  plate  ( N  m'2  ).  The 
tangential  stresses  on  the  top  of  the  acoustic  coating  are  modeled  as  free  boundary  conditions  and 
are  written  as 


rzy(x,y,c,t)  =  p  2 


dvj  (x,  y,  c,  t)  dw2  ( x,y,c,t ) 


8z 


dy 


=  0 


and 


Tzx(x,y,c,t)  =  p2 


au2(x,.y,c,/)  | 

dx  dz 


=  0. 


(24) 


(25) 


The  acoustic  pressure  in  the  fluid  medium  is  governed  by  the  three-dimensional  wave 
equation  and  is  written  in  Cartesian  coordinates  as 


d2pa(x,y,z,t)  |  g2p„(x,y,z,/)  a2p„(x,y,z,Q 
ax2  ay2  az2 


1  d1  pa(x,y,z,t) 


'f 


dt‘ 


=  o. 


(26) 


where  pa(x,y,z,t)  is  the  pressure  (  N  m~“  )  and  c f  is  the  real-valued  compressional  wavespeed 

in  the  fluid  (ms'1).  The  interface  between  the  fluid  and  the  top  of  the  acoustic  coating  (z  =  c ) 
satisfies  the  linear  momentum  equation,  which  relates  the  normal  acceleration  of  the  plate 
surface  to  the  spatial  gradient  of  the  pressure  field  by 


d*-w2(x,y,c,t)  a  pa(x,y,c,t) 

p<  - = - Tz - 


where  p  y  is  the  density  of  the  fluid  ( kg  m'3 ) 


(27) 
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3.  ANALYTICAL  SOLUTION 


All  three  displacements  are  now  determined  by  solving  equations  (8),  (9),  (10),  and  (11), 
inserting  these  results  into  equations  (5),  (6),  and  (7),  and  applying  equation  (4)  to  eliminate  two 
of  the  resulting  constants.  Because  of  the  periodicity  in  the  x-direction,  the  functional  form  of 
the  displacements  are  equal  to  the  sum  of  an  unknown  function  in  the  z-direction,  which  is 
multiplied  by  an  indexed  exponential  function  in  the  x-direction,  by  an  exponential  function  in 
the  y-direction,  and  by  an  exponential  function  in  time.  Others  note  this  functional  form  and 
discuss  its  derivation.6  This  is  written  as 

777  =  +  00 

u j(x,y,z,t)=  Y  U w  * (z) exp(i£wx) exp(i/r vy) exp(-i att ) ,  (28) 

777  = -00 


777  —  +0C 

vj(x,y,z,t)=  Y  ^i/\z)exp(i£wx)exp(i£vy)exp(-i^),  (29) 

m——cc 

and 

777  =  4-00 

wj(x,y,z,t)=  Y  ^m  \z)^p(ikmx)exp(ikyy)exp(-\a)t)  (30) 

777 — —oo 


where  i  =  V-T ,  co  is  frequency  ( rad  s'1  ),j  is  either  1 ,  which  corresponds  to  the  backing  plate  or 
2,  which  corresponds  to  the  acoustic  coating,  kv  is  wavenumber  with  respect  to  the  y-axis 

( rad  m’1 ),  and 


km=kx  + 


2mn 


(31) 


where  kx  is  wavenumber  with  respect  to  the  x-axis  ( rad  m'1 ).  The  functions  u\j\z) ,  vjrf\z) , 
and  W^j)(z)  are  given  by 
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(32) 


um\z)  =  Am){km  cosla^z]  +  B^\km  sin  [a^z] 


+ 


C(/) 


m  (7ysin[/?iy)z]  -  D</}  ^cos[/#>z]  + 
Pm  Pm 

k; 


pU) 

^  m 


?U) 


/*; 


0) 

w 


sin[^)z]-F^) 


A 


U>  +  • 

w  ~(y) 


A 


m 


cos[^J/)z] , 


vm  ' )(z>  =  4/^v  cos[arJ/)z]  +  F^iA:  sin[ctr^/)z]  + 


CU) 


m 


nU) 

Mm  ^  /  ;  x 

n(j) 
Mm 


sin[/^/)z]  +  D 


(./') 

/?? 


A 


(./) 


v  w 


m 


Mm 


cos  [y?J/)z]  + 


_  /r(./) 


n  +  Fj/>  cosl^'z] , 

Pm 


m  / 
n(j) 
Mm 


(33) 


and 

^'(z)  =  )a«)  sin[Qri/’)z]  +  5^')a^/)  cos[«J„7)z]  + 

-  r'/’iA-  ,.  cos^z]  -  D^'xky  sin [/?</>z]  +  (34) 

+  Em]'lkm  cost/^z]  +  FY,)ikm  sin[^/}z] . 

where  ,  B^\  C^;),  D^/\  and  F^7)  are  unknown  complex  wave  propagation 
coefficients  of  the  backing  plate  and  acoustic  coating,  a\p  is  the  modified  wavenumber 

( rad  m  1 )  associated  with  the  dilatational  wave  and  is  expressed  as 

am]  =^[kdj)]2-kl  ~-k2y  ,  (35) 

where  A:^7  *  is  the  dilatational  wavenumber  on  the  layer  (/')  and  is  equal  to  col ;  j3Yp  is  the 
modified  wavenumber  ( rad  m  * )  associated  with  the  shear  wave  and  is  expressed  as 

AJ)  =  >]2  “*v  .  (36) 

where  kYl)  is  the  shear  wavenumber  on  the  layer  (/')  and  is  equal  to  col . 
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The  pressure  field  in  the  fluid  consists  of  an  incoming  incident  pressure  wave  applied  to  the 
structure  at  the  top  (z  =  c)  and  an  outgoing  scattered  wave  that  propagates  in  the  positive  z- 
direction.  Using  the  analytical  form  for  the  pressure  field  of 

/w  =  +oo 

p(x,y,z,t)  =  (z)  exp(i£wx)  exp(\kyy)  exp(-i  cot ) ,  (37) 

m=- oo 


results  in  the  solution  of  equation  (26)  where 

pm(z)  =  prn  e*P(irmz)  +  (>moPl  exp(-i ymz),  (38) 

where  SmQ  is  the  Kronecker  delta  function,  Pj  is  the  magnitude  of  the  applied  pressure  field 
(Nm‘  ),  and  ym  is  the  modified  wavenumber  ( rad  m  )  associated  with  the  compressional 
wave  in  the  fluid  and  is  given  by 


The  relationship  between  arrival  angles  on  the  plate  of  an  acoustic  wave  and  the  x-  and  y- 
wavenumbers  is  determined  by 

kx  =  (<y/cy  )sin(#) 
and 

kv  =(a>/c  /  )sin(^) . 

where  6  is  the  arrival  angle  of  the  acoustic  wave  (rad)  with  respect  to  the  x-axis  and  (/>  is  the 
arrival  angle  with  respect  to  they-axis  (rad).  A  value  of  0=  <f>=  0  corresponds  to  a  broadside  or 
boresight  wave  exciting  the  system.  For  the  remainder  of  the  paper,  the  exponential  function 
with  respect  to  time  is  surpressed  in  the  equations. 


(40) 


(41) 
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The  displacement  fields  in  equations  (28)-(30)  are  now  inserted  into  the  boundary  value 
equations  given  by  equations  ( 14)  through  (25).  Additionally,  the  pressure  field  in  equation  (37) 
and  the  interface  equation  listed  as  equation  (27)  are  utilized,  and  the  resulting  equations  become 


/W  =  +00 


4  X  iM4!)(«)exP(i*m*)  +  'lli*y  X  ^/ml)(«)exp(i^m.v) 


+ 


n=+  oo 


m=+  oo 


2]  ^,i)(a)exp(i^OTx) 


S(x-nL ),  (42) 


m  -  +oc 


A77  =  +O0 


/q  X  — ^7 — exp(i£wx)  + /q  £  ff//!,l)(a)exp(i^wx)  =  0 , 


m  =  -oo 


m  =  -cc 


(43) 


W  =  +00 


W  =  +00  f  (1)  /  X 

/*!  X  i^w^,>(«)exp(i^wx)  +  ^1  ^ - -exp(iAwx)  =  0 , 

m=-cc  m=—oo 


dz 


(44) 


_  AZ7  =  -t-OU 

4  X  ^m^m'C^expO^  +  A,^  £  rj,1)((>)exp(i£m.v)  + 


m=+  oo 


m=+oo 


(A, +2//,)  ]T  ~~  exp(i^wx’)  =  A2  £  i^mt/i2)(6)exp(i^wx)  + 


dz 


m=+oc 


m=+oo 


X  ^2)(^)exp(i^m.Y)  +  (A2 +2//2)  Y,  ^T-^exp(i^x), 


(45) 


w=+oo  3 ^(1)  ( b )  m-+oo 

M\  X  — — exp(i£mx)  +  /qi£v  £  ^(/^exp^x)  = 


w=+oo  QV^\b)  m=+ 00 

A2  X  — ^7 — exp(i*mx)  +  //2i*v  £  fFj,2)(6)exp(i£mx), 


(46) 
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m=+ oo 


M\  X  i*m^«)(fc)exp(i*w*)  +  /*1  X  (/))exp(i£wx)  = 


m-~oc 


dz 


rfl  —  +QC  m  =  -»-Q 

/*2  X  ^m^m2)(^)exp(i/:OTx)  +  z/2  X 


m=~oo 


'^expOW, 

<9z 


w  =  +oo 


w  =  +oo 


X  t/i!)(^)exp(i^)=  2  £/i,2)(6)exp(i£wx), 


m=-cc 


m=- oo 


/?7  =  +00 


w=+oo 


X  ^l)(/))exp(i^wjr)=  51  ^(^expU^x), 


m=-oo 


m=- oo 


/w=-i-qo 


W7  =  +00 


X  ^1)(/>)exp(i^x)=  X  ^2)(6)exp(i^wx), 


m=-oo 


m=-  oo 


m=+oo 


^2  X  i^mt/'w2)(c)exp(iA:wx)  +  /l2i^y  X  Pi2)(c)exp(i*,„x)  + 


m=+cc 


0*2+  2/y2  )  X  exP(^mx )  +  X 


m=+a>f <y2/?/  ^ 


^  y 


^2)(c)exp(i^mx) 


=  -2P/  exp^A^x) , 


W  = 


/*2  X 


^dV^ic) 


m=+  oo 


m=- oo 


_w _ 

dz 


-exp(i/rwx)  +  //2^v  X  ^i,2)(c)exp(Uwx)  =  0, 


m=-oo 


and 


W  =  +00 


m  =  +CO  fir y (2)  /  \ 

M2  Z  ^w^2)(c)exp(i/:wx)  +  /y2  £  — ®- - exp(i^x)  =  0, 

m=-oo  m=-oo  UZ 


(47) 


(48) 


(49) 


(50) 


(51) 


(52) 


(53) 


13 


To  eliminate  the  Dirac  delta  comb  function  present  in  equation  (42),  a  form  of  the 
Poisson’s  summation  formula  written  as 


«=+c o  j  n=+  oo 

^S(x-nL)-—  ^exp(i2;zraVL) 

n  =  —oo  n = —oo 


(54) 


is  used.  Furthermore,  the  identity 


■  n=+  oo 

-  X 

L  ^ 

n  =  -oc 


m  =  +  oo 

Z  ^,!)I)(«)exp(iA'wx) 

exp(i2<wv/L)  =  — 

n  —  -Too 

I 

_TYl  —  — OO 

L 

_fl — — oo 

m  =  +co 


£exp(i*wx)  (55) 


m  =  ~oo 


is  also  applied  to  equation  (42).  This  modified  version  of  equation  (42),  with  equations  (43)- 
(53)  are  all  multiplied  by  exp(-i£px)  and  integrated  from  [0,Z.].  Because  the  exponential 
functions  are  orthogonal  on  this  interval,  the  equations  decouple  into  sets  of  w-indexed  equations 
expressed  as 


X\\kmU^(a)  +  A\i kyV^l\a)  +  (T]  +  2//] ) 


(1). 


dlV^ia) 


K_ 

L 


n=+  oo 

fl  =  —00 


(56) 


dV^ia)  .,  m 

»\—*L  l  +  M\'kyW%\a)  =  Q, 

cz 


(57) 


”>(a)=o, 

CZ 


(58) 


i, !*„,[/<]>(*)  +  2,i* ,,V$\b)  +  (2,  +  2/i, )aW”)W  - 

«/i2)  (*)  +  ^i*  v421  (»)  +  U2  +  2//2 )  , 


(59) 
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(60) 


P\ 


Ml WS' (*)  +  Al -f~-  =  //2 i*m^2) (ft)  +  A2  <?C/”'(/’>  . 

dz  dz 


U{Jl\b)  =  u£\b), 
v£\b)  =  v£\b), 
w£\b)  =  w£\b). 


(61) 

(62) 

(63) 

(64) 


^2  xkmuffl  (c)  +  A2  i  k  y  V i2)  (c)  +  (i2+2/i2)  CW^)(C) 

<7Z 


+ 


C02pf ' 

wl?\c)  =  . 

ir"  J 

-2  P, 

0 


m  =  0 
m*  0 , 


(65) 


^2 


+  Pi  ikyWffl (c)  =  0 , 


and 

//2itm^2»(c)  +  A2^|^  =  0. 

cz 


(66) 


(67) 


The  functional  form  of  the  displacements  given  in  equations  (32),  (33),  and  (34)  are 
inserted  into  equations  (56)-(67)  and  the  resulting  algebraic  matrix  equation  for  each  w-indexed 
coefficient  is 
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[A(Aw)]{x(Aw)}  =  2jF(A-„)]{x(*„)M 


0  m  *  0  , 


(68) 


where  [A(AW )]  is  a  twelve  by  twelve  matrix  that  models  the  dynamics  of  the  backing  plate  and 


the  acoustic  coating,  {\{km )}  is  the  twelve  by  one  vector  of  unknown  wave  propagation 
coefficients,  [F(A„)]  is  the  twelve  by  twelve  matrix  that  models  the  dynamic  interaction  of  the 
springs  and  the  backing  plate,  and  Jp)  is  the  twelve  by  one  matrix  that  models  the  incoming 
acoustic  wave  acting  on  the  structure.  The  entries  of  the  matrices  and  vectors  in  equation  (68) 
are  listed  in  the  appendix.  Equation  (68)  is  now  written  for  all  values  of  the  index  m  and  the 
results  are  inserted  are  rewritten  in  global  matrix  form  as 

A  x  =  F  x  +  p,  (69) 

where  A  is  a  block  diagonal  matrix  and  is  written  as 


A  = 


A(A_,)  0  0 

0  A(*0)  0 

0  0  A  (A,) 


(70) 


F  is  a  rank  deficient,  block  partitioned  matrix  and  is  equal  to 


F(A_])  F(A0)  F(A] ) 


F  =  F(*_!)  F(*0)  F(A, ) 

F(A_,)  F(A0)  F(A|) 


(71) 
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p  is  the  system  excitation  vector  and  is  written  as 


P  = 


(72) 


and  x  is  the  vector  that  contains  the  unknown  wave  propagation  coefficients  and  is  equal  to 


x  = 


(*<*-, )>T  W*o)lT  M*i)!t 


where 


{x(^o)} = 


j(l)  d(0  /^(l) 

Ao  °o  Lo 


D0)  £(»  F0) 

u0  c0  r0 


.A 


(2)  b(2)  r(2)  n(2)  r(2) 


0 


0 


c, 


o 


D, 


0 


'0 


A\) 


(73) 


(74) 


The  0  term  in  equation  (70)  is  a  twelve  by  twelve  matrix  with  all  zero  entries  and  the  0  term  in 
equation  (72)  is  a  twelve  by  one  vector  with  all  zero  entries.  The  solution  to  the  wave 
propagation  coefficients  is  now  found  by  truncating  the  matrix  equation  (69)  to  a  finite  number 
of  terms  and  solving 


x  =  [A-F]-'p.  (75) 

Once  these  are  known,  the  displacement  field  of  the  system  in  the  spatial  domain  can  be 
determined  using  equations  (32),  (33),  and  (34).  Furthermore,  the  stress  distribution  in  the 
backing  plate  and  the  scattered  acoustic  field  can  be  computed. 
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4.  MODEL  VALIDATION 


The  fully  elastic  model  that  has  been  developed  in  sections  2  and  3  can  be  compared  to  a 
fluid-loaded,  ribbed,  Bemoulli-Eular  thin  plate  that  has  been  previously  developed.7  4  This  will 
provide  validation  of  the  model  for  low  frequencies  and  small  plate  thickness.  This  Bemoulli- 
Eular  model,  however,  only  incorporates  flexural  wave  behavior,  making  the  model  assumptions 
invalid  at  higher  frequencies  and  wavenumbers.  The  thin  plate  model  has  one  degree-of-freedom 
that  is  constant  displacement  in  the  z-direction.  This  displacement  equation  is  written  as 


m=+  oo 


w(x,y)=  £fTmexp(i£wx)exp(i/:v>;)  , 


(76) 


and  this  displacement  field  can  be  determined  by 


(»}=[T+Rrljf), 


(77) 


where 


M  =  {-  W_x  fV0 


(78) 


T’(-l)  o  0 


[T]  =  •••  o  r(0)  o 
o  o  r(l) 


(79) 
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1 


[R]  =  (A7L) 


and 


{f}=  {•••  0  -2  Pi  0 


In  Eq.  (79),  the  indexed  entry  is 


T{n)  =  D(k4n+k2nk;+k*)-phco2+(pfco2  n r„). 


(80) 


(81) 


(82) 


where 


D  = 


E7?3 

12(1  -o2)  ' 


(83) 


where  E  is  Young’s  modulus  (  N  m  )  and  v  is  Poisson’s  ratio  (dimensionless). 

Figure  2  is  a  plot  of  the  magnitude  of  the  transfer  function  of  the  plate  normal  displacement 
( w )  divided  by  the  amplitude  of  the  applied  incident  pressure  field  ( Pj )  versus  k  v  wavenumber 
(top)  and  ky  wavenumber  (bottom)  at  a  frequency  of  10  Hz.  In  the  top  plot,  the  kv 

wavenumber  was  equal  to  2  rad  m  *  and  in  the  bottom  plot  the  k  x  wavenumber  was  also  equal 

to  2  rad  m'1 .  This  extremely  low  frequency  was  chosen  because  it  is  a  value  at  which  the  thick 
plate  model  should  theoretically  agree  with  the  thin  plate  model.  This  example  was  generated 

using  the  following  system  parameters:  total  plate  thickness  h  is  0.005  m,  plate  density  p  is  1200 
kg  m'3 ,  Lame  constants  and  X2  are  9.3 1  x  1 08  N  m'2 ,  Lame  constants  px  and  p2  are 
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8  2  3 

1.03x10  N  m'" ,  fluid  density  pj  is  1025  kgm"  ,  fluid  compressional  wavespeed  c^is  1500 
m  s'1 ,  rib  stiffness  per  unit  length  K  is  2  x  104  N  m'“ ,  spatial  location  x  is  0,  spatial  location  v 
is  0,  and  rib  separation  distance  L  is  0.5  m.  For  the  thick  plate  model,  the  plate  interface  region 
location  b  is  -0.003  m  and  the  output  location  was  z  =  -h/ 2,  which  is  -0.0025  m.  In  figure  2,  the 
solid  line  is  the  fully  elastic  plate  theory  developed  above  and  corresponds  to  equation  (30)  and 
the  dot  symbols  is  the  Bemoulli-Eular  plate  theory  and  corresponds  to  equation  (76).  The  elastic 
plate  model  was  truncated  to  seven  modes  that  produced  an  84-by-84  system  matrix  and  the  thin 
plate  model  was  truncated  to  5 1  modes.  There  is  almost  total  agreement  between  the  two  models 
over  both  wavenumber  variables. 
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Figure  2.  Transfer  Function  Magnitude  of  Normal  Displacement  Divided  by  Incident 
Pressure  at  10  Hz  with  ky  =  2  rad  m'  (top)  and  kx  =  2  rad  m~'  (bottom)  for  Thick  Plate  System 

( _ )  and  Thin  Plate  System  (  •) 
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5.  AN  EXAMPLE  PROBLEM 


An  example  problem  is  now  generated  and  discussed.  This  problem  consists  of  a  0.0065-m 
(0.25-inch)  thick  aluminum  plate  coated  with  a  0.0254-m  (1-inch)  thick  urethane  polymer 
coating.  This  example  was  generated  with  the  following  parameters:  aluminum  Lame  constant 
Tj  is  5.1 1  x  1010  N  m'"  ,  Lame  constant  p\  is  2.63  x  10*°  N  m’2 ,  density  pj  is  2710  kg  m’"*, 
urethane  Lame  constant  is  2.09  x  109  N  m'2 ,  Lame  constant  is  7.14  x  107  N  m"2 , 
density  pi  is  1 1 10  kg  m'  ,  urethane  structural  loss  factor  is  0.05,  fluid-compressional 
wavespeed  cy  is  1475  ms'1 ,  fluid  density  pf  is  1025  kg m'3,  internal  rib  spacing  L  is  1.00  m, 
and  rib  stiffness  per  unit  length  K  is  1 .20  x  1010  N  m'2 .  The  problem  is  investigated  from  two 
technical  standpoints:  (1)  the  reception  of  an  acoustic  signal  and  (2)  the  echo  reduction  of  an 
acoustic  emission  directed  at  the  system.  Because  the  major  interest  of  this  analysis  is  response 
to  an  acoustic  wave  in  the  fluid,  subsequent  figures  are  shown  with  respect  to  arrival  angles 
rather  than  wavenumbers.  The  relationships  between  arrival  angles  and  wavenumbers  are  given 
by  equations  (40)  and  (41). 

The  reception  of  an  acoustic  signal  in  the  urethane  can  be  accomplished  by  an  internal 
array  of  sensors  embedded  in  the  coating.  A  typical  sensor  in  a  solid  material  will  detect  average 
normal  stress,  determined  mathematically  using  the  equation 

_  ,  x  rxv(x,y,z5,ry)  +  f,3,(x,y,z5,ry)  +  f2Z(x,y,z5,ry) 

rave(jc,y,z5,<y)  = - - - - - - - ,  (84) 

where  zs  is  the  location  of  the  sensors  along  the  z-axis  (m).  Figure  3  is  a  plot  of  the  magnitude 
of  the  average  stress  divided  by  the  incident  pressure  versus  #  arrival  angle  for  ^arrival  angle 
equal  to  0  degrees  (top)  and  average  stress  divided  by  the  incident  pressure  versus  #  arrival  angle 
for  ^arrival  angle  equal  to  45  degrees  (bottom)  at  a  frequency  of  8027  Hz,  a  value  ofy  =  0,  and  a 
value  of  z.v  =  -0.0102  m.  In  figure  3,  the  solid  line  is  the  fully  elastic  plate  theory  with  ribs  at  x  = 
0,  the  long  dashed  line  is  the  fully  elastic  plate  theory  with  ribs  at  x  =  LI 2,  and  the  short  dashed 
line  is  the  fully  elastic  plate  theory  without  ribs  at  any  jc  value.  Note  that  the  addition  of  the  ribs 
has  made  the  stress  field  become  spatially  varying,  and  that  different  magnitudes  will  be  present 
at  different  values  of  jc. 
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Figure  3.  Transfer  Function  Magnitude  of  Average  Stress  Divided  by  Incident  Pressure  at 

8027  Hz  with  <f>  =  O°(top)  and  <f>  =  45  ° (bottom)  for  the  System  with  Ribs  atx  =  0  ( _ )  ,  x  = 

L/2  ( _ )  and  System  without  Ribs  at  any  x  Value  ( - ) 
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Once  the  average  stress  field  has  been  calculated,  the  summed  acoustic  response  to  an 
incident  pressure  field  at  fixed  frequency  and  (f>  =  0  (or  k  v  =0)  can  be  computed  using 

N 

B(kx,(o)  =  ^rave(x„ ,0,z5,ft>)exp(-iA:5xw) ,  (85) 

n= 1 

where  N  is  the  number  of  sensors  in  the  array,  ks  is  the  steered  wavenumber  of  the  array  with 
respect  to  the  x-axis  ( rad  m"1 ),  and  xn  is  the  location  of  the  «th  sensor  (m).  The  effect  of  the  ribs 
is  examined  here  with  respect  to  a  16-element  linear  array  oriented  in  the  x-direction  with  an 
interelement  spacing  of  0.0625  m.  Using  these  parameters,  the  summed  acoustic  response  of  the 
array  is  computed  using  a  uniform  shading  function  and  the  results  are  shown  in  figure  4.  In 
figure  4,  the  top  plot  is  the  array  steered  to  0°  and  the  bottom  plot  is  the  array  steered  to  -45°. 

The  solid  line  is  the  array  response  for  the  system  with  ribs  and  the  dashed  line  is  the  array 
response  for  the  system  without  ribs.  Note  that  for  the  array  steered  to  0°  and  -45°,  the  rib 
structure  slightly  widens  the  beam  and  surpresses  the  maximum  response.  Additionally,  the  side 
lobe  structure  of  the  array  response  is  shifted  outward  to  a  greater  angle  for  the  ribbed  response. 
The  ribs  also  surpress  the  energy  in  the  side  lobes  past  the  second  side  lobe. 
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Figure  4.  Summed  Acoustic  Array  Response  with  a  Steer  Angle  of  0°(top)  and  -45° (bottom) 
for  the  System  with  Ribs  ( _ )  and  System  without  Ribs  ( - )  at  802  7  Hz 
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In  underwater  structures,  the  scattered  pressure  field  is  an  important  quantity  for  quiet 
operation.  The  scattered  pressure  field  divided  by  the  incident  pressure  field  is  given  by  the 
expression 


P  (  X  V  0)\  Wl  —  +<X)  p  , 

~  p  : —  =  exp(i^x)exp(i^y)+  £  - — ^^2)(0)exp(i^wx)exp(i/:v>;)  ,  (86) 

where  Ps(x,y)  is  the  scattered  pressure  field  in  the  spatial-frequency  domain,  the  first  term  on 

the  right-hand  side  is  the  reflected  pressure  field,  and  the  second  term  on  the  right-hand  side  is 
the  radiated  pressure  field  caused  by  normal  displacement  of  the  plate.  Figure  5  is  a  plot  of  the 
magnitude  of  the  scattered  pressure  divided  by  the  incident  pressure  versus  <9  arrival  angle  (top) 
and  <f>  arrival  angle  (bottom)  at  a  frequency  of  8027  Hz.  In  the  top  plot,  the  </>  arrival  angle  was 
equal  to  0°  and  in  the  bottom  plot,  the  ^arrival  angle  was  equal  to  45°.  In  figure  5,  the  solid  line 
is  the  fully  elastic  plate  theory  with  ribs  at  x  -  0,  the  long  dashed  line  is  the  fully  elastic  plate 
theory  with  ribs  at*  =  LI 2,  and  the  short  dashed  line  is  the  fully  elastic  plate  theory  without  ribs 
at  any  *  value.  The  ribs  add  about  a  2  dB  spatial  variation  to  the  scattered  pressure  field  when 
compared  to  the  system  without  ribs.  This  differential  is  frequency  dependent  as  it  is  more 
pronounced  at  various  frequencies. 
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Figure  5.  Transfer  Function  Magnitude  of  Scattered  Pressure  Divided  by  Incident  Pressure 

at  802  7  Hz  with  <f>  =  O°(top)  and  <j>  =  45°  (bottom)  for  the  System  with  Ribs  at  x  -  Q  ( _ )  ,  .v 

=  L/2  ( _ )  and  System  without  Ribs  at  any  x  Value  ( - ) 


28 


6.  CONCLUSIONS 


A  fully  elastic  analytical  model  of  a  system  that  consists  of  a  ribbed  plate  covered  by  a  fluid- 
loaded  acoustical  coating  has  been  derived.  The  model  has  been  shown  to  agree  with  previously 
developed  thin  plate  solutions  to  the  problem  at  low  frequency.  An  example  problem  was 
developed  to  illustrate  high-frequency  behavior  and  a  comparison  between  the  ribbed  and 
unribbed  structures  is  included.  The  specific  solutions  of  embedded  sonar  system  performance 
and  scattered  pressure  field  were  investigated.  It  was  shown  that  the  ribs  have  a  moderate  effect 
on  the  high  frequency  dynamic  response  of  the  structure. 
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APPENDIX 

MATRIX  AND  VECTOR  ENTRIES 


The  entries  of  the  matrixes  and  vectors  in  equation  (75)  are  listed  below.  Without  loss  of 
generality,  the  top  of  the  top  plate  is  defined  as  z  =  c  =  0.  For  the  [A  (k„ )]  matrix,  the  nonzero 
entries  are 


al,l  =  WltC^j1*)2  +lin  +^y)]_2//l(clr«1))2}cos(««1)a)’ 

(A- 1 ) 

a\,2  =  {-*l[(an ))2  +kn  +^v)]-2/y1(ai1))2}sin(or^1)a), 

(A-2) 

«1,3  =2Wl^1)^sin(^1)a), 

(A-3) 

«1,4  =  — 2Wl  Pn)ky  cos(fin^a)  > 

(A-4) 

a\,5  =  ~2'M\  Pn)kn  sin  , 

(A-5) 

"1,6  =  ^M\Pn)kn  COS {p^a)  , 

(A-6) 

a2,\  =-2\fi\a^)ky  sinCar^a), 

(A-7) 

a2,2  =2iV\an)kycos(ay)a), 

(A-8) 

°2,3=-M(til))2  +kn  -k2y]cos(/J^a), 

(A-9) 

*2,4  =-/^i[(^1))2  +  *2  -A^sinO^a), 

(A- 10) 

A-l 

"2,5  =-2p\kykncos(/3^)a),  (A- 1 1 ) 

a2,6  =~2V\kykn  sin(/3(nl)a),  (A-12) 

"3,1  =-2ifiJ\al})ktlsin(al1])a),  (A-13) 

a32  =2i{/\ain])kncos(aj1])a),  (A-14) 

"3,3  =^l^«^v[cos(/?i1)a)-sin(^1)a)],  (A-15) 

"3.4  =M«Msin(^il)")-cos(/9<1)")],  (A-16) 

"3.5  =M\  [(Am')2  -  k},  ]cos(^,1)a)  -  /J]kj,  sin(/3^a) ,  (A-17) 

"3.6  =Mi [(fin]))2  -k*]sin(filil)a)-finky  cos^^a),  (A-18) 

"4,1  =  {-^[("i0)2  +kl  +^v)]_2y"l(an  ))2}cos(«i1)6),  (A-19) 

"4.2  =  WitC^i0)2  +kn  +A:v)]-2//1(aJ1))2}sin(aJ1)Z)),  (A-20) 

"4.3  =2iM\/3n)kySm(^)b),  (A-21) 

"4.4  =  — 2i//i P^ky  cos (/3{„])b),  (A-22) 


A-2 


*4,5  =-2i^l  Pn)kn  sin  (fiPb), 

*4,6  =2[/J\Pn)kn  COS(/3{„])  b) , 

*4,7  =  U2[(«i2))2  +  ^  +^)]  +  2//2(a®)2 

*4,8  =^2[(«i2))2  +  +kj)]  +  2M2(aP)2 

*4,9  =  -2'/^2  tf sin(^2)/>)  , 

*4,10  =  2i//2  fin2)ky  cos(p(2)b), 

*4,11  =2\n2p£)k„sm  (p(2)b), 

*4,12  = -2i>*2^2)^  cos  (fii2)b)t 
*5,1  =  -2i,*i««  v  sinCa^/?) , 

*5,2  =  2i*i«J%-  cosCorJ^)  , 

*5.3  =  ~M\[(Pl!))2  +  *2 -k^cos^b), 
*5,4  =-/^l[(^i1))2  +  *2 -k2y]sin(ril)b), 
*5,5  =  ~2JU\kykn  cos(fi{n])b) , 


(A-23) 

(A-24) 

}cos  (a(2)b). 

(A-25) 

}sin  (a{2)b). 

(A-26) 

(A-27) 

(A-28) 

(A-29) 

(A-30) 

(A-31) 

(A-32) 

(A-33) 

(A-34) 

(A-35) 

A-3 

(A-36) 


"5,6  =  -2//i  kyk„  sin (P^b), 

"5,7  =  2'M2ah2)ky  sin(or,(2)/>) ,  ( A-37) 

"5,8  = -2\<U2ah2) k y  cos(a{„2) b) ,  (A-38) 

"5,9  =  M2[{P\2))2  +kl  ~k;.]cos(p{2)b) ,  (A-39) 

"5,10  =M2 [(fih2))2  +k%  -k2y]s\n(P(2)b),  (A-40) 

"5,11  =2  V2kyk„cos(/3(2)b),  (A-41 ) 

"5,12  =  2 M2kykn  sin (P(2)b) ,  (A-42) 

"6,1  =  -2i//1aJ1)A:wsin(a(,1)6),  (A-43) 

"6,2  =2i//1or,()1)^cos(«,(,l)/7),  (A-44) 

"6,3  =  V\k„ky[cos(/3{„])b)-sin(pl,l)b)],  (A-45) 

"6,4  =^i^«^[sin(^1)6)-cos(/?i1)6)],  (A-46) 

"6,5  =  V\[(Pn])2  -k2n]^pPb)-nxklsm{P^b),  (A-47) 


A-4 


*6,6  =M\[(Pn))2  -  k2  ]sin(/?^6)  -  ju\k2  cos(fi^b) , 

(A-48) 

*6,7  =  2i /J2a(2)kn  sin(«J,2)/)) , 

(A-49) 

*6,8  =  -2i^2«i2)^«  cos(or  (,2)/>) , 

(A-50) 

*6,9  =  -Miknky[cos{/l(2)b)  ~  s\n(^2)b)] , 

(A-5 1 ) 

*6,10  =  -y*2*«Msin(/^2)6)  -  COS(/?J2)/>)] , 

(A-52) 

*6,11  =-A2[(^2))2  -^]cos(^2)6)  +  //2^2sin(^2)6), 

(A-53) 

*6,12  =->*2  [(^2))2  -  k2  ]sin(^2>Z>)  +  /i2ky  cos  (/}{2)b). 

(A-54) 

*7,1  =i^„cos(aJ1)6), 

(A-55) 

*7,2  =  •*«  sinCaJ^ft). 

(A-56) 

*7,3  =  (1f  srn^^fc), 

(A-57) 

"7'4=  in  cos(aV ’*)• 

dS0 

(A-5  8) 

°7.5=  A',' '  +  n,  sin 

(A-59) 

A-5 


"7,6  ~ 


"7,7  = 


"7,8  - 


"7,9  = 


"7,10  “ 


"7,1 1  = 


"7,12  = 


"8,1  =  ' 


"8,2  =  * 


"8,3  - 


k 2  ^ 

Pm 


cos  (fi^b), 


-i kn  cos (a[^b) , 


-\kn  sin(«j2^) , 


-k  k 

n  v  •  //?(2)x,\ 

— — -sin(/?;  ]b ), 

A."* 


I  b 

KnKv (2), 


Pyn 


(2) 


•cos(^'A), 


fiKn 


(2) 


*v  ^ 

A',2' 


;,2  > 

^2)  +  V 


(2) 


cos(/?,(;2)/?) , 


ky  cos (a^b) , 


A:v  sin(a),1*/?), 


.2  3 

o(l)  _^7_ 

Pn  /?(1)  , 
k’n  J 


sinip^b). 


(A-60) 

(A-61) 

(A-62) 

(A-63) 

(A-64) 

(A-65) 

(A-66) 

(A-67) 

(A-68) 

(A-69) 


A-6 


( 


cos  (P^b), 


(A-70) 


"8,4 


"8,5 


"8,6 


fiS'  + 


.2  } 


Pi" 


"  ‘  sin {p^b). 


ft 


(1) 


Mv 


ft 


(i) 


cos  (ftl)b). 


(A-71) 


(A-72) 


"8,7  =-ikyCos(a{„2)b), 

(A-73) 

"8,8  =  -i^v  sin(a^2)6), 

(A-74) 

(  k2  \ 

"8,9=  /^)+  — 737  sin(/?i2)6), 

l  ti2)) 

(A-75) 

f  k2  \ 

"8,10=  P?  +  /L  COS(/?^)6), 

l  ’  J 

(A-76) 

«8.1l  =  ,2|  SITU/?,1,-'*), 

«2> 

(A-77) 

°8,12  =  COS  (0(„2>b), 

tf1 

(A-78) 

A-7 


"9,1  =  sin (a^b) ,  (A-79) 

"9,2  =  a?' 1  cos(«J1)/>) ,  (A-80) 

"9,3  --iky  cosijSj,l)b),  (A-81) 

"9,4  =-i^sin(^1)A),  (A-82) 

"9,5  =  cos(/?^6) ,  (A-83) 

a96=\knsm(/3^)b),  (A-84) 

"9,7  =  «»2)  sin(aj2)6) ,  ( A-85) 

"9,8  =  ~««2)  cos(arJ2)f>) ,  (A-86) 

"9.9  =  i£  v  cos (J3(„2)b) ,  (A-87) 

"9,10  =i^vsin(^2)6),  (A-88) 

"9,11  =~il<ncos(/3(n2)b),  (A-89) 

"9,12  =  ->*#»  sin (P(n2)b) ,  (A-90) 

"10,7  =-h[(ai,2))2  +  *«  +^v)]-2/y2(«J2))2,  (A-91) 


A-8 


(A-92) 


"10,8  = 


(0 


PfCC 


(2) 

n 


[r„ 


"10,9  = 


-arpfk^, 

Yn 


"10,10  ="2i  HlP^ky, 


co  pfk„ 

"10,11  = - . 

Yn 

"10,12  =  2iM2  fil2)kn, 

"n,8  =2i  p2ai2)ky, 

"11,9  =  ~M2 (fih2)  +kn  -ky), 
"11,11  =-2y"2  knkyi 

"12,8  =2i/i2«i2)^, 

"12,9  =M2knky> 


"12,10  =  -M2knky, 

01241 


(A-93) 

(A-94) 

(A-95) 

(A-96) 

(A-97) 

(A-98) 

(A-99) 

(A- 100) 

(A- 101) 

(A- 102) 

(A- 103) 

A-9 


and 


"12,12  =~M2  ky- 

The  nonzero  entries  of  the  [F(£„)]  matrix  are 
/l,l  =~(K/L)  a{nl)  sinCa^o) , 

/l,2  =(K/L)  a\ 7n  cos (a(„])a) , 

f\,3  =~(K/L )  \ky  cos (p^]a) , 

/l,4  =  ~(K/L)  iky  si n(j3{n])a), 

f\,5  =(K/L)ik„  cos {/3{„])a), 
and 

/l,6  =(K/L)ikn  sin (/?Jna). 

The  entries  of  the  p  vector  are 

p  =  [0  00000000-  2P, 


(A- 104) 

(A- 105) 

(A-106) 

(A- 1 07) 

(A-108) 

(A-109) 

(A-l  10) 

0  0]T.  (A-l  1 1 ) 
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